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Abstract 



The explicit formulae for m-soliton solutions of (l+2)-dimensi- 
onal matrix Davey-Stewartson equation are represented. They are 
found by means of known general solution of the matrix Toda chain 
with the fixed ends [0 . These solutions are expressed trough m + m 
independent solutions of a pair of linear Shrodinger equations with 
Hermitian potentials. 
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1 Davey— Stewart son Equation 

Let u, v are two non-singular s x s matrix functions of x, y, i.e. each matrix 
element is a function of x, y coordinates of two-dimensional space. Partial 
derivatives of these functions up to some sufficient large order are assumed 
to exist. 

We define the matrix Davey-Stewartson equation (DSE) as the following 
partial differential equation. 

!Kl + ^ + ^- 2oM /^- 26 /^-)„.^0 (1 ) 

where a, b are arbitrary real numbers and z* is a Hermitian conjugate of a 
matrix z. It will be convenient to deal not with the equation (TO), but with 
the following expanded system, which we call the matrix Davey-Stewartson 
system (DSS). 



iut + au xx + buyy — 2au J dy{vu) x — 2b J dx(uv) y • u — 

(2) 

— ivt + av xx + bvyy — 2a J dy{vu) x ■ v — 2bv J dx(uv) y = 

Bellow for definity we had chosen a = b = 1. It is easy to see that the DSE 
is the system (|2|) under additional condition v = u*. This condition we call 
the condition of reality. 
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In the case s — 1 (scalar case), when u,v are scalar functions and the 
order of the multipliers is not essential (|l|) is the usual, well-known Davey- 
Stewartson equation 0. In the scalar case soliton solutions of the DSS were 
obtained in ||. 



2 Discrete substitution 

The method we use to solve the problem is based on the discrete transforma- 
tion investigation. Here we consider concrete discrete transformation, which 
is important for our problem. 

By direct calculations can be checked that is invariant with respect 
to the following change of unknown matrices u, v. 



u — v , v = vu — (v x v ) y v = v uv — {v 1 v y ) x (3) 

Here u,v denote the "new", transformed operators. Invariance means that 
matrices u, v satisfy exactly the same system which matrices u, v satisfy. 
Mapping (|3|) is an invertible one and the "old" matrices u, v can be expressed 
trough the "new" ones. 



fi- 1 , 



u 



UV — {UyU V ) 



u = u 



VU - (U l U X )y 



(4) 



Transformation (H) can be rewritten in the form of an infinite chain of equa- 
tions in two equivalent ways as 



U 



n+1 



(5) 



or as 



V n -l v n ~ V n l V n+l 



(6) 



where (v n , u n ) is a result of n-times substitution @ applied to some initial 
matrices v ,u . Sequences (^]|J) with vz\ = vn = boundary conditions we 
call the matrix Toda chain with the fixed ends. 

In the scalar case s = 1 general solution of the Toda chain with the fixed 
ends was found in M for all series of semi-simple algebras, except E7,E$. In 
[||] this result was reproduced in terms of invariant root technique applicalible 
to all semi-simple series. 
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The explicit general solution of the matrix Toda chain with the fixed 
ends was found in . It was expressed trough N + N arbitrary independent 
matrix functions of a single argument X r (x), Y r {y) as 

N 

V = J2 X rYr (7) 
r=l 

To (0) correspond the following formula for u^. 

u N = jr Yr (x)X r (y) (8) 

r=l 

Matrices X,Y here are not arbitrary ones, but in some way depend from 
X, Y. Both these results ([?]) and (§) will be used in further consideration. 

3 General Strategic 

We are going to solve the DSS (|2|) under condition of reality u — v*. Here we 
describe how the discrete transformation is used for that. General idea is the 
following. At first we take some obvious solution of the DSS (Q). It may be 
not a solution of the problem (reality condition may be not satisfied) . Then, 
by means of the discrete transformation (^), we get from that initial, obvious 
solution a solution which is satisfy to the condition of reality. 

For uo = the first equation of the system (|2|) is satisfied identialy, the 
second one gives 

- iv 0t + v 0xx + v 0yy + Vi(t, x)v + v V 2 (t, y) = (9) 

where V\ , V2 are arbitrary s x s matrix functions of their arguments (these 
terms arise from the undefined integrals / dx(uv) y , J dy(uv) x in the system 
@)- Obviously, condition of reality is not satisfied for this solution. But 
after enough times discrete transformations ([3]) it is possible to come to the 
solution for which it is satisfied. To clarify this, let us consider some solution 
u, v for which condition of reality is satisfied, u = v*. Denoting u\,vi and 
U-i,V-i the results of direct @ and inverse @ substitutions respectively, it 
can be easily proved that U-i = V\* and v-i = u\. On the m-th step, we 
have: w_ m = v m * and V- m = u* m) where index m (— m) correspond to the 
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result of the m-times direct (inverse) transformation. And vice versa one can 
prove that if we begin from solution uq = 0,vq and after 2m-times discrete 
transformation receive u 2m = ^o,^2m = 0, the solution in the middle of the 
chain automatically satisfy the reality condition, u m+ \ = v^ +1 . 

The system arising from equations u = u 2m = is already resolved by 
the formula (^). So it remains to solve the equation: 'u 2m = v o- ^ leads to 
the following relations between X r ,X a M and Y r ,Y a [ r j 

X* = X CT [ r ] Y* = Y a [ r ] (10) 

where a denotes one of the (2m)! possible permutations of the 2m low indexes. 
To solve (|K]) at first it is necessary to find the dependence of X and Y from X 
and Y respectively. Finally equation (Q) in terms of X r , Y r can be rewritten 

as 

- iX rt + X rxx + Vi(t, x)X r = - iY rt + Y rxx + Y r V 2 (t, y)=0 (11) 

Thus to find the m-soliton solutions of the DSE ([]]) it is necessary to under- 
take the following steps 

• find the dependence Xj (Xl, . . . , X 2m ) and % (Yi, . . . Y2 m ) 

• solve the system (|I0~D 

• find such a dependence of matrix functions X r , Y r from the time argu- 
ment, which will satisfy to the system (|TT|). 

After this substituting X r ,Y r in (^) we find v , for which u m+ i = v^ +1 is 
some partial (m-soliton) solution of the Davey-Stewartson equation ([!]). 



4 Scalar Case 

To gain some experience firstly we consider the scalar case s = 1, for which 
much of the necessary calculation steps are well-known and much simpler 
then in the general matrix case. 

In this case for the mentioned above boundary conditions the following 
formulae for arbitrary k takes place [0]. 

Detk-i Detk+i 
u k = n . v k = n j. > Det-i = 0, Det = 1 (12) 
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where Detk is the principle minor of dimension k of the matrix (v° = Vq) 



(v 



° v l 

y xy 

yy %yy 



XX 

V° 

^ xxy 

J° 

xxyy 



\ 



V 



/ 



and v° is determined by (|7|) , where X r , Y r are arbitrary scalar functions of 
their arguments. Substituting (|7|) in the formula for w 2m from (|12|) and 
comparing with (H) we find 

W / 2m-l(^l) X2, X r -i, X r+ i, ...X 2m ) 



(13) 



W2m(X\ 1 X<2i ....X2m) 

Here and bellow Wk denotes a Wrosnkian of dimension k constructed from 
the functions in the brackets. 



Wk(0i, ■■■,9k) 



9i 
9[ 



92 
92 



(fc-1) (fc-1) 
9i 92 



9k 
9l 



9 { r ] 



W 



(14) 



Expressions for Y r can be received from (|13"D by the simple exchange X — » Y. 

In the condition of reality ( |lOl) we use the permutation cr[r] = 2m — r + 1 
To resolve ([[(]) and (|TT|) it is suitable to represent the functions X r , Y r in the 
Frobeniouse-like form 



From 



we find 

2m 



Xi = 4>i, X r = (fiij dx(f> 2 
Y 1 = i/j 1 , Y r = Jdxip2 ■ 

2m 



dx(f) r 



(15) 



X 2m = (II (j>k) X r = (II ^fc) / dx<p2m ■■■ f 

k=l k=l J J 

After this the reality condition flTU| ) takes the following form 



dxcj) 



2m~r 



■ > 2m—r-\- 



2 (r = 2, 3,. ..2m) 



J m+1 



J m+1 



n mi 

k=l 



(16) 



(17) 
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And from (|TT| ) we have 



r-l 



k=l 



(18) 



The imaginary unity i here is included into the time variable, which therefore 
should be treated as a pure imaginary from this moment. Independently can 
be checked that the systems (|I7|) and ( jT8|) are compatible and if (|18|) is obeyed 
for some <p r , r < m, for 2m -r+2 it also holds. Therefore it is enough to 
consider only equations with r < m in the system (0). Now we introduce 
the new unknown functions j~ x — - ■ ■ 4> r , r < m + 1. From (|18D we find 



(/-7r-l)t = " (/- 1 /r-l(ln/ r /,- 1 )') / 



(19) 



From (0) it follows that = / m+1 . Substituting this in the (m + l)-th 
equation of the last system, we have 



(fmfn 



fmf m (l n /m/ m _i 



(20) 



Equation (|20|) is equivalent to one-dimensional Shrodinger equation with 
arbitrary real potential. 



f mt + C = Uf m , u = u* 
Now let us consider the m-th equation of the system (|T9|). 

{fmfm-l)t = — (/ m 1 /m-l(ln/ m / m -l)' 

Partially resolving it as 

fmfm-l = Z', (lnf m f,, 



Zt 



I m—1 1 



(21) 



(22) 



(23) 



and excluding function / m _i, we conclude that function zf m satisfy exactly 
the same equation ([21]), which f m satisfy. Denoting u i: (1 < i < m) m 
independent solutions of (|2~T|) we find 



fm ^1) zfn 



U 2 



m—1 



2 /rr 



■Ul ^2 



"1 



(24) 
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In the general case, for arbitrary % the following formula holds. 

fr = r < m (25) 

m—r 

where Wi = Wi(u\, . . . ,Ui). 

To prove ( p^D we use the well-known Jacobi identity for determinants. Let 
T is some infinite in both directions matrix; D n (T) denotes the determinant 
of its n x n principle minor, T s is the matrix received from T by deleting 
its s-th column and T p -by deleting its £>-th row. In this notations Jacobi 
identity takes the form 

D n (T)D n (T") - D n (T n )D n (T n ) = D n+1 (T)D n ^(T) (26) 

From (OH) the following identity can be easily derived. 

WiW' t - WiWt = W^W i+1 (27) 

where Wi = Wi(ui -> u i+l ) = W^Ut, u 2 , . . . , Mj_i, 

Now let us partially resolve equation (|T9| ) for arbitrary r. 

fr-l = z'f r = - Z j (28) 

Excluding f r _i from the last system we find that f r and zf r are different 
solutions of the same equation. And if f r is given by the formula (p5|) , zf r 

can be determined as 

W m - r +i 



Zfr 



W m _ r 



After this from (p8|) with the help of the identity ( p7|) we easily find 

"m-f+2 



/r-1 



m— r+1 



Thus formula ( p5|) is proved by induction. 

Finally for functions <p r from the definition of f r and formulae flTT] ) we 
have 

, * , Wm-r^Wm-r 

(p m+ l = VlV{ (p r = 7772 

' * m— r+1 

(29) 

W m -\ 

01 = -TTr <K = 4>2m-r+2 T < 771 

Win 



The analogues expressions take place for functions ipk 



W m - r+ 2W m - T 

(30) 



w 2 . , 

" ' m—r+1 



4>1 = —TT r iPr = V2m-r+2 T < Tfl 

In ( PU| ) Wi = Wi(vi, . . . , Vi) and u 4 = Vi(y) 1 < i < m are m independent 
solutions of (1 + 1) -dimensional linear Shrodinger equation with some arbit- 
rary real potential V. 

v it + v? = Vvi V — V* (31) 



5 Matrix case 

Here we consider a general problem as it was formulated in the sections 1 
and 3. We find m-soliton solutions of the DSE for an arbitrary dimension 
of the unknown matrix u. We therefore receive matrix generalizations of all 
formulae of the previous section. It turns out that quasi-determinant s of 
matrices with non-commutative entries play the role of usual determinants. 
Conception of a quasi-determinant was introduced recently by Gelfand and 
Retarh 0. We use an independent technique, more appropriate for our 
particular case, but quasi-determinants can be used as well. 

With the chain (|5]||) under mentioned above boundary conditions we 
connect the following recurrent relations 

n-l 

R n = v^v ny , fl£ = £ (Si; 1 + R.Sf 1 ) (32) 

fc=0 

with the boundary conditions Sf = 1 for arbitrary i. From definitions (|32l) 
and equations (HH) we easily find 

n 

Sl n = H R n S = V^Voy ...y (33) 

9 

v n+1 = -v n (Sl +1 ) x = (-lr+'voiSlUSl), ■ ■ ■ (Sl +1 ) x (34) 
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For such introduced ( |32"D matrix functions S% the following relation is true. 



qq 



(35) 



Now let us find the dependence of X from X. For that we use the fact 
that each matrix function X$ is determined only by matrices Xi, . . . ,X 2m 
and therefore we can choose matrices Y\ , . . . , Y 2m in an arbitrary way. It is 
convenient to choose 



i-l 



Y; 



+ 



y 



2m- 1 



(2m - 1)! 



X 



2rn 



(36) 



where E is a unity sx s matrix. Substituting ( 53 ) in the expression for Vi m -\ 
from (03) we find 



7 2m-l 



y=0 



^(^(^•••(r^ 



where matrices T" are determined by the following relations 



(37) 



rpq 



with the boundary conditions 



(TLl)x _1 (Tn-l) 



rpq _ Q q 

1 o — °0 



y=0 



(38) 



(39) 



Expression (^) correspond to one of the functions Xj. And because these 
functions can be enumerated in various ways we can choose 



X\(Tq) x (TI) x ■ ■ ■ {T 2 \ 



2m-2)x 



F(xi, ...,X 2r 



(40) 



Formula for arbitrary i can be received from ( [40]) by the cycler permutation 
of the indexes. 



X 



-i 



F(a i [x 1 , . . .,x 2m ]) 



(41) 



Arbitrary multiplier can be added in the formula (0). It will be counted 
in the expression for Yi. We added (— 1)* _1 to do further calculations more 
convenient. After this using fl4~T| ) and (|3"1"D we find 



Yi 1 = ~(Ql 



2m-2)y 



(Ql) y Y 1 = G(Y 1 ,...,Y, 



2m 



(42) 
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Yr^ = (-iyG(a l [Y 1 ,...,Y 2m ]) 



(43) 



where 



\{Qn-X 



gs = n t . 1 >r 1 ( 44 ) 

Now as in the previous section we represent the initial functions X r ,Y r in 
the Frobenious-like form. 



(45) 



X 1 = 4> 1 , X 2 = fa Jdxfa, X 3 = fafdxfafdxfo, ... 
Y 1 = v/j 1 , Y 2 = Jdxfa ■ fa, Y 3 = Jdxfjdxfa-fa^j -fa, 



After permutation of the indexes the formulae for X r will coincide with ([lq) . 
The only difference is that in the matrix case the order of the multipliers 
must be taken into account. 



Xt=p, X 2 = dx(j)2m -P, X 3 = / dx / dx(f)2m-l • 02m " P, 



Yi = -S, Y 2 = -Sj dxip2m, Y 3 = -Sj dxip 2m J dxip 2m -l, ■ ■ ■ 

(46) 

where p = (fa • • ■ 4> 2m ) 1 and s = ("02m • ' • fa) 1 ■ The condition of reality 
taken in the form X r = X*, Y r = Y* gives 

<K = fam-r+2 2 < T < 171, 
4>m+l = (0m+l) _1 = (0102 



J rn+l) — lViV2 ' ' ' Vm) ' " " ^mj 

0* = ^2m-r+2 2 < T < T71, 

V'm+l = {fam+lfa 1 = ' ' ' V>l)*(VW0m-l ■■■fa) 



(47) 



The fact that all functions X r are solutions of the same equation (|il[ ) leads 
to the following system. 

- (fa) t + (2(0 1 02-0 s -i)~ 1 (0i02-0 s -i) / s + = (48) 
Introducing new the functions j~ x = fafa ■ ■ ■ 4>r from (H8|) we find 



(fr-lf r X )t — fr-lfr 1 — fr-l(f r ^ 



(49) 
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Then from fl4"T|) and ( pj[ ) we conclude that / m is a solution of (1 + 1)- 
dimensional linear Shrodinger equation with Hermitian potential. 

fm t + f^ = Wf m , W = W* (50) 

Solution of the system ( JI5| ) can be found by the same scheme as in the 
previous section. Matrix case does not require to use the Jacobi identity, 
because recurrent definitions are used. 

h = u x f m - r = (U 1 r )'---(Ul)'u l (51) 

Matrix functions U% are determined by the following recurrent relations. 

u q n = {u£\)' [(ul,)']' 1 

with the boundary conditions 

where matrices u r are different solutions of the equation fl5H|). From this we 
find the formulae for (p m - r for arbitrary r. Finally we have 

01 = A _1 4> m - r = {U})' 0<r<m-2 

(52) 

(j) m+ l = U-iU* (j)* r = (p2m-r+2 2 < T < 171 

ip m - r = {U})' < r < m - 2 (53) 
ip m+ i = -v x v\ ip* = i) 2m -r+2 2 < r < m 

V3 



For -0j we find 



where 



and matrices Vi(y) are different solutions of (1 + l)-dimensional linear 
Shrodinger equation with arbitrary Hermitian potential. 

v it + v'l = Vi M, M = M* 

Now substituting ([521, 53) directly in (|45j) and (0) and then in the formula for 
v m+ i from ([34]) we find the m-soliton solution of the matrix DSE. We do not 
write down the corresponding expression, because it can be easily received, 
but is too large to represent it here. 
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6 The Simplest Example of One— Soliton 
Solution 

Substituting m — 1 in the formulae of the last section we find 
v = X X Y X + X 2 Y 2 , X 1 = (j> 1 , X 2 = (j) 1 j dx<p 2 

Yi = ipi, Y 2 = J dxifj 2 ■ ipi 

After this we find the following expression for the one-soliton solution of the 
matrix DSE. 

Mi = Vf 1 (l + f dx(f> 2 J dyip^J fa 1 (54) 

Matrix functions <f>(t,x),if)(t,y) are determined by u,v solutions of the one- 
dimensional linear Schrodinger equations. 

01 = M^ 1 , (j) 2 = UU* 

ipi = v' 1 , ip 2 = —v*v 

u t + u xx + uMi(t,x) = 

vt + v yy + M 2 (t, y)v = 0, Mi, 2 = M* 2 



7 Conclusion 

The main result of the paper is the explicit expressions for the m-soliton 
solutions of the (l+2)-dimensional matrix Davey-Stewartson equation. By 
means of the corresponding formulae of the sections 4 and 5 these solutions 
are expressed trough m + m independent solutions of a pair of linear (1+1)- 
dimensional Shrodinger equations. 

From the group-theoretical point of view it means that we had realized 
the finite-dimensional representation of the group of integrable mappings. 
This viewpoint remained beyond our concrete calculations. 

Note that restriction with the finite-dimensional matrices is absolutely 
nonessential. We had never used this restriction and moreover the dimension 
(s) was not included in any expression. 
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